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This work is a sequel to our work "The Spin Density Matrix I: General Theory and Exact Master Equations" 
(eprint cond-mat/0708.0644). Here we compare pure- and pseudo-spin dynamics using as an example a system 
of two quantum dots, a pair of localized conduction-band electrons in an n-doped GaAs semiconductor. Pure- 
spin dynamics is obtained by tracing out the orbital degrees of freedom, whereas pseudo-spin dynamics retains 
(as is conventional) an implicit coordinate dependence. We show that magnetic field inhomogeneity and spin- 
orbit interaction result in a non-unitary evolution in pure-spin dynamics, whereas these interactions contribute 
to the effective pseudo-spin Hamiltonian via terms that are asymmetric in spin permutations, in particular, the 
Dzyaloshinskii-Moriya (DM) spin-orbit interaction. We numerically investigate the non-unitary effects in the 
dynamics of the triplet states population, purity, and Lamb energy shift, as a function of interdot distance and 
magnetic field difference AB. The spin-orbit interaction is found to produce effects of roughly four orders of 
magnitude smaller than those due to AB in the pure-spin model. We estimate the spin-orbit interaction mag- 
nitude in the DM-interaction term. Our estimate gives a smaller value than that recently obtained by Kavokin 
[Phys. Rev. B 64, 075305 (2001)], who did not include double occupancy effects. We show that a necessary and 
sufficient condition for obtaining a universal set of quantum logic gates, involving only two spins, in both pure- 
and pseudo-spin models is that the magnetic field inhomogeneity AB and the Heisenberg interaction are both 
non-vanishing. We also briefly analyze pure-spin dynamics in the electron on liquid helium system recently 
proposed by Lyon [Phys. Rev. A 74, 052338 (2006)]. 



PACS numbers: 

I. INTRODUCTION 

The spin degree of freedom of a locaUzed particle, e.g., an 
electron or nucleus, is a popular carrier of quantum informa- 
tion. It serves as a qubit which can be manipulated in order to 
accomplish a computational task. The spin of electrons local- 
ized in quantum dots (QDs) or by donor atoms has been the 
subject of extensive recent studies ' .2-3-4i5,6,7,8,9, lo, 1 1 , 12, 13, 14. 15 ^ 

Consider two electrons trapped in two sites A and B, e.g., 
two QDs each containing one electron. The two-electron sys- 
tem is fully described by the total wavefunction |5'tot), which 
depends on the electrons' coordinates r and spin variables cr. 
The two-electron spin-density matrix, obtained by tracing out 
the orbital degrees of freedom, p = Tr^? j'l'tot) (4'tot|, fully 
describes the spin dynamics as long as one cannot or does not 
wish to apply measurements that can separate or localize elec- 
trons spatially; the only observable is then the electron spin, 
s" = ^ctq, where ctq are the Pauli spin one-half matrices with 
a — x,y, z. Since the spin system is not closed - there is a 
coupling to the electrons' spatial degrees of freedom - we ob- 
serve open system effects, i.e., the spin dynamics becomes in 
general non-unitary. We refer to this dynamics as pMre-spin 
dynamics. 

In contrast, pseMiio-spin dynamics is the standard case 
where the electron spin observable is not free from coordi- 
nate dependence but includes information about the electron's 
localization orbital. In the /jieMc/o-spin case one defines the 
electron spin operator as a bilinear combination of electron 
annihilation and creation Fermi operators, cas, c^^., in a lo- 
calized orbital (1>a (s is a spin index, A is the QD index): 



= 5 I]L'=i ^As i'^a)ss' CAs', a = x,y, z. Then the oper- 
ators {sJIIq obey the usual su(2) commutation rules. 

This paper is the sequel to our work Ref. [l^ (henceforth 
"part I"), where we derived an operator-sum representation 
(OSR) as well as a master equation in the Lindblad and time- 
convolutionless (TCL) forms for the spin-density matrix of a 
two-electron system. In this sequel we focus on a detailed 
comparison of pure and pseudo-spin dynamics. We are in- 
terested in particular in how non-unitary effects in pure-spin 
dynamics are translated into the corresponding unitary ones in 
pseudo-spin dynamics and vice versa. We show that as long 
as there is no magnetic field inhomogeneity the pure-spin dy- 
namics is unitary, but in the presence of magnetic field inho- 
mogeneity this dynamics is non-unitary 

The paper is organized as follows. We begin, in Section 
|ll]by highlighting the differences and relationship between 
pseudo- and pure-spin models. Sectionlllllprovides a concrete 
illustration in terms of a system of two QDs trapping one elec- 
tron each. In it, we examine the role of different interactions 
in both pseudo- and pure-spin dynamics. We first derive the 
coordinate part of the Hamiltonian (subsection IIII Al l and the 
form of the dipolar interaction (subsection IIII Bl i. In subsec- 
tions lIII Cl and llll Dl respectively, we then present calculations 
illustrating effects due to both external magnetic field inhomo- 
geneity and the spin-orbit interaction in the pure-spin model. 
In subsections IIII El and IIII Fl we discuss universal quantum 
gates in both pseudo- and pure-spin models. Subsection llllGI 
presents our estimates for spin-orbit interaction effects in the 
pseudo-spin model, and compares these estimates to the re- 
sults of Ref. 17. We conclude in SectionlVl 
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Atomic units, h = e = rUe = 1, 1/c ~ 1/137, are used 
throughout the paper unless stated otherwise. 



II. PSEUDO- VS PURE-SPm APPROACHES 

In this section we discuss the relation between the present 
approach based on the spin-density matrix and the pseudo- 
spin effective Hamiltonian approach. The latter is usually de- 
veloped as a low-energy mapping within the Hubbard model 
Hamiltonian of interacting electrons^' " ' ^ '9,20,2 1 .22,23 ^ 
not follow the Hubbard model since it is highly simplified and 
neglects many interactions which we would like to keep here. 
For Hubbard model analyses in the quantum computation con- 
text, see, e.g., Ref.[^ 



A. Pseudo-spin effective Hamiltonian 

In order to keep the present treatment as simple as possi- 
ble we restrict ourselves to the two orbitals approximation 
used in part I; inclusion of excited-state orbitals is straight- 
forward. Consider the four single-occupancy basis states 
{<i><;i, $ti, i = 1,2,3}, where $si is a singlet wavefunction 
with two electrons localized on different QDs, A and B, while 
^ti are the corresponding triplet wavefunctions. The two 
double-occupancy states {$s2- ^ss} describe two electrons in 
a singlet state, localized on the same QD, A or B. The total 
wavefunction (t) in this basis set takes the form 



where 



n,s{E) = PHP + PHQ 



1 



E-QHQ 



QHP. 



(3) 



Observe that Eq. (|2]i is exact but non-linear, and has 6 solu- 
tions. 

Due to interelectron repulsion the double occupancy states 
are usually much more energetic than the singly-occupied 
ones if the electrons are well localized in QDs. We consider 
the low-energy physics described by Eq. (|2]l where the total 
energy E is near the energies of singly-occupied states. In 
general Hcs is not a Hamiltonian since it is a function of the 
energy E. However, if the energy gap between the P- and 
Q-states is large enough, one can expand and approximate 



(E — E] 

n,s{E) = Hcs{E) + y PHQ—^ L 

t^L {E - QHQ) 

/(I) 



-QHP 



n,ffiE)+nll>iE) {E-E)+ 0[{E . 
n,s{E)+n^^^{E) {E-E), 



E) 



(4) 



where 

/(I) 



E is an average energy in 

Hl'f^iE) = PHQ {E - QHQ)~^ QHP 



the P-subspace and 
Keeping terms up 



to the first order in Eq. dUi, the non-linear Eq. (|2| can be 
reduced to a generalized Unear equation problem 



1=1 

where the complex amplitudes {asi{t), ati{t)} define, respec- 
tively, the singlet and triplet states population. In the total 
Hilbert space, the state is defined by 11 real parameters [12 
real parameters defining {asi{t),ati{t)} minus a normaliza- 
tion condition]. The unitary evolution in the total Hilbert 
space is described by 



\at{t)) 



exp {—iHt) 



Wsio)) 



where H is the total two-electron system Hamiltonian. 

Since these basis states are orthonormal, projection opera- 
tors into the corresponding subspaces can be written as 



Q = |$,,2) ($.2! 
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El 



s3| 



(1) 



where Q projects onto the double occupancy states. Then, 
using the method of projection operators, one obtains the 
Schrodinger (eigenvalue) equation projected into the P- 
subspace 



(nes{E)+n'^^J{E)E - (1 + nl}J{E))E'^ = O. (5) 

Solving Eq. (|5]), we obtain four low energy solutions; the two 
high energy, double-occupancy solutions are lost in this ap- 
proximation. Therefore, in the low energy, pseudo-spin ap- 
proximation the state is described by 7 real parameters. 

In the following, we assume n'-^^ (P) EE for simplicity. 
The effective Hamiltonian Eq. (01 can be recast into a pseudo- 
spin form. Using Eq. ([T]) we have 



3 

+Y (^f l^si) {<^u\ + nr m , (6) 



where 



(*,l|7^eff(P)|*.l) 



W** = (<ft.|7^off(P)|$t,), 



(7) 



(7^cff(P)-P)P* = 0, 



(2) 



In the second quantization representation, the P-subspace ba- 
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sis vectors take the form 



1 / 



V2 I'^T'' 

^(IT>A®U>B-U>A®lt>B), 



„t 



|<fti) - c\^cl^\Q) = \])^®\])s, 



1 



^(IT>A®U)i3 + U)A®lt>B). 
1**3) = ct^^4jO) = |i)^®U)5, 



(8) 



where we introduced pseudo-spm states |s)^, s =ti = 
A, B, localized near the A and S sites [the term pseudo em- 
phasizes the fact that these are not really spin states since 
they depend on the electron orbital degrees of freedom]. 
Eqs. ([8]l establish a one-to-one correspondence between 4 ba- 
sis states {$si , ^ti, i — 1,2,3} and 4 tensor-product pseudo- 
spin states \s)^(^\s')p, where s, s' —1, |, a, /3 = A, B. Then, 
relabeling the ;?seM(io-spin states as |0, 1) = |t, i) and intro- 
ducing the /?ieMt/o-spin Pauli and identity operators 

a, = |0)(1| + |1)(0|, 

ay = -z(|0)(l|-|l)(0|), 

a. = |0)(0|-|1>(1|, 

/ = |0)(0| + |1)(1| (9) 

where we temporarily dropped the subscripts A and B, one 
easily finds that 



= S, |<fH)($t,| 



= K,. (10) 
Here the psewJo-spin operators S, Tij are defined as 

S = jI-sa-sb, T^ii ^ \I + \Sy^ + sazSbz, 

T22 = -^I + SAxSBx + SAySBy — SAzSBzi 
T33 — jl — ^Sz + SAzSBz, 

T12 = [^S+ + Js] , T23 = ^ [^5*+ - Js] , 

Ti3 = SA'xSBx — SAySBy + i (sAxSBy + SBxSAy) , 

-.t 



T3I = Tig, 

where 



L32 



^23 



(11) 



Js = SAzSBx + SAxSBz + i{sAzSBy + SAySBz) 

S± = Sx±iSy, S^sa + sb, 
and K is defined as 



Ki = 



2V2 



K 



3 — 



^ /— 1 1 Jas ) ^ M Ja. 

2V2 \ yy 



(12) 



where Jas ~ sb — sa >^ S . In fact, Eqs. ( fTTT i and ( fT2] i can 

be obtained from the corresponding ones in part I if the pure- 
spin operators si_2 are replaced respectively by the pseudo- 
ones, sa.b- We reproduce these formulas here in order to 
make the presentation as self-contained as possible. 

As is seen from Eqs. (fTTT i and (fT2l i. the first line of Eq. (|6]l is 
symmetric with respect to spin permutations [A <-> B], while 
the second one is asymmetric representing, in particular, the 
Dzyaloshinskii-Moriya (DM-type) interaction termi2i2^ No- 
tice that these asymmetric (in spin permutations) terms cancel 
out of unitary spin dynamics after averaging over orbital de- 
grees of freedom, as demonstrated in part I. However, they do 
not disappear completely, but rather are converted into the cor- 
responding non-unitary terms plus the Lamb shift term as will 
be seen in next subsection. From the symmetric part of the 
Hamiltonian (|6]l, using Eqs. (fTTT i one can derive the isotropic 
Heisenberg exchange interaction term 



Ti-H — Jrsa ■ SB, 



(13) 



where the Heisenberg exchange interaction constant Jh — 
i Hll—H^"; in contrast, as was demonstrated in parti, the 
Heisenberg interaction term does not affect the unitary evolu- 
tion of the spin density matrix, apart from the Lamb-energy 
shift. In subsection nil CI we demonstrate numerically the ef- 
fects of the Heisenberg interaction on both the Lamb-energy 
shift and the non-unitary part of the spin density matrix evo- 
lution. 

Observe that the asymmetric part of the Hamiltonian Eq. 
(|6]l is proportional to the singlet-triplet subspace interaction 
matrix Hf, which is responsible for the coupling between sin- 
glet and triplet states. As will be demonstrated in Section HID 
the non-zero coupling between these states is due to B-field 
spatial inhomogeneity (i.e., it cannot arise due to the homo- 
geneous component of the external magnetic field), as well as 
due to the spin-orbit interaction. 



B. Spin density matrix 

In part I we derived the Lindblad-type master equation for 
the spin density matrix 



dpjt) 
dt 



(14) 



ij ^ ' iJ 



H" -J P 



where the first and second terms describe, respectively, uni- 
tary and non-unitary contributions to the evolution. H" is 
an effective pure-spin Hamiltonian which includes the Lamb- 
shift term, jPq; the pure-spin operators and K,; are de- 
fined by Eqs. (fTTT i and ( fTST i where sa.b ^ Si 2- The index 
a — {s,t,m} specifies which initial state p(0), singlet, s, 
triplet, t, or a mixed one, m, is taken. 
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As mentioned in part I, all the matrix functions in Eq. (fl4] i 
as well as the pseudo-spin Hamiltonian Eq. (|6]l are expressible 
in terms of H matrix elements 



H = 



a, P = s, t, 



(15) 
i,j = 1,2,3. 



In the following example, we consider the triplet case for 
which we have 



xl = i{Qo 

Pa. — ~^ 

where 



(16) 



fe 



Csfc) {esk\Ra{0) + \esk) {etkW x 
X;exp(-i£fci) [\etk) {esk\ i?a(0) + \etk) {etk\] ] 

k J 

(17) 

Here, Ra{Q) is a correlation matrix, which establishes an ini- 
tial correlation between the singlet and triplet amplitudes 



a,(0) =i?„(0)at(0) 



(18) 



[in the triplet case, we have i?Q=t(0) = 0; in the mixed case, 
where both a,(0) ^ and Of (0) ^ 0, i?a=m(0) ^ 0; for the 
singlet case, see part I] and e^, \esk), and \etk) are solutions 
to the eigenvalue problem 



^SS JJSt 
JjtS Jjtt 



\esk) 
\etk) 



£k 



\esk) 
\etk) 



■ ,6. 
(19) 



III. EXAMPLE: SYSTEM OF TWO QUANTUM DOTS 

In this section, we investigate the role of different inter- 
actions in the calculation of the H matrix. Let us consider 
a system of two electrons trapped at sites ta and re (rA,B 
are radius-vectors of the centers of QDs in the z = plane) 
created by a system of charged electrodes in a semiconduc- 
tor heterostructure so that the electrons are confined in the 
z = plane or a system of localized conduction-band elec- 
trons in n-doped GaAs as in our calculation example. The 
heterostructure trapping potential 



(20) 



is separable in the in-plane and out-of-plane directions; V± (z) 
and VA.B{r) are the trapping potentials in the z-direction and 
in the z = plane around fA,B respectively. If the elec- 
tron system is placed in a constant magnetic field Bq directed 



along the z-axis (with vector potential 



fx Bq 



), 



then the in-plane motion, in a superposition of the in-plane 



confining oscillatory potential and a perpendicular magnetic 
field, is described by the Fock-Darwin (FD) states.— Approx- 
imating the confining potential by a quadratic one 



(21) 



we can take as basis "atomic" orbitals the ground-state func- 
tions 



<l>A,B{z,'r) = (po{z)RA^B{\f- fA,B\) 



(22) 



where the out-of-plane motion in the z-direction is "frozen" in 
the ground state (Pq{z) in the potential V±{z), and the ground 
FD state is 



R 



FD 

A,B 



1 



^ttIa.b 

Ir 



exp 



4/2 

^'■a.b 



(23) 



1 



1^1 



Here Ia.b is the effective length scale, equal to the magnetic 
length Ic in the absence of the confining potential, uja,b = 0; 
ijJc = Bq/c\s the cyclotron frequency. 

The orbitals Eq. ( l22b must be orthogonalized. One way 
to do this is a simple Gram-Schmidt orthogonalization proce- 
dure: 



DA = 

Ib = 



{4>B - S4>a) , 



where the overlap matrix element Sab 
{(f'A \4'b) can be calculated analytically 

„ '^IaIb 



S 



BA 



72 I 72 

I A + l-B 



exp 



AB 



(24) 



S 



(25) 



For appropriate values of system parameters such as the in- 
terdot distance tab and the external magnetic field Bq, the 
overlap becomes exponentially small. 

The other, more symmetric way is to make a transition to 
the "molecular" or two-centered orbitals by pre-diagonalizing 
the coordinate part of Pauli's non-relativistic Hamiltonian he, 
which describes the electron's motion in a superposition of the 
trapping potential and magnetic fields: 



4>B 



CAA(PA + CAB(PB, 
CBA't'A + Cbb4>B; 



A,B, 



(26) 



The two-state eigenvalue problem Eq. ( |26l ) is solved analyt- 
ically in terms of "atomic" orbitals matrix elements: hij — 

In general, given the "molecular" Eq. ( |26] l or "half- 
molecular" Eq. ( l24b basis choices, one cannot ascribe a spin 
to a particular QD, since an electron in a molecular orbital 
belongs to both QDs. 

The total Hamiltonian contains both coordinate and spin- 
dependent terms. First we consider the coordinate part of the 
Hamiltonian in the (I)a,b basis set. 



5 



A. Coordinate part of the Hamiltonian 

In view of the orthogonality of the singlet and triplet spin 
wave functions, the spin-independent part of the Hamiltonian 
does not contribute to the singlet-triplet coupling, = 
Hl^ = 0, whereas for the singlet-singlet and triplet-triplet 
Hamiltonians we get 



we have -ffjjp 



TTts 

-"dip 



and a non-zero contribution to 



i/^fp comes only from the contact term: 



(i/|fp)., = (1.45-2^)meV(/..|,5(ri2)|/s,) 

dii di2 di3 
1.09AmeV | ^22 



il3 2 



-^33 



(31) 



rjss 

-«cll - 
rrss 



H. 



cl3 



ttss _ 
-«c22 - 



TTSS 

-«c33 

St 

where 



f TTSS \3 TJSS^ 

\ cij S i^j — l") cij cji 

hAA + hsB + Vcc{AB; AB) + Sco(AB; BA), 
V2 (Iba + iUAB; AA)) , (27) 

V2 {jiAB + «oc(AB; BB) 

2hAA + VcciAA; AA), H'J^ = icc{AA; BB), 
2hBB+iee{BB;BB), 



£tl, 

hAA 



JiBB +Vcc{AB;AB) 



(28) 

VeeiAB;BA), 



hij = {(t>i h, 



i,i,k,l = A,B 



I '^c '■Fj / 7 ^7 j A, B 
1 



£^■12 



Mmm), (29) 



with hij = Eidij for "molecular" orbitals and -Doc being the 
interelectron electrostatic interaction matrix elements. The 
matrix H^'^ is diagonally dominated if the overlap S* ^ 1; 
H^li is the singlet energy of the singly occupied state whereas 
H^22 ^iid -ffc33 ^re energies of doubly occupied states if one 
neglects the coupling between single- and double-occupancy 
states. Observe that the Heisenberg constant Jh = St — £s 
where is the lowest eigenvalue of the matrix H^'' . The ma- 
trix elements hij and Vcdij', kl) can trivially be expressed in 
terms of the corresponding matrix elements h^j and Woc(u; kl) 
where the orthonormalized states c^/s are replaced by 0/s us- 
ing the relations Eq. (|24] | or ( |26] |. 



where A is an effective constant of the interaction that confines 
electrons in the z-plane and 



dii 
di3 

d22 



1 

;2 I ;2 "^^P 
'■A + 'S 



' AB 



V2(4>A 1 4) 



^2 0/2 , ;2 

/^U/Ib 
V 2 777T5 — — exp 



4 3?|+/i 



' AB 



4 311 + 1% 



4iy 



-'33 



4/2 ■ 



(32) 



The magnetic dipole contribution to the triplet-triplet interac- 
tion Hamiltonian can be written as 



i/^*p = 0.36 meV 



/ to 

3 



3 J* 

775^1 -2to :7^tl 



V 



-3*2 -^h 



to 



(33) 



where 



'2 71 

0, 1, 2 are dipole tensor operators 

1 - 3 C0S2 012 



to 



tl = 



to 



sin 26*12 exp(i(/?i2) 



(34) 



' 12 



sin^ 9i2 exp(2i(y9i2) 



' 12 



with (ri2, 6*12 , 'P12) being spherical coordinates of the inter- 
electron radius-vector fi2 — ri — r2', the bar over ti denotes 
averaging over the triplet coordinate wavefunction: 



t^ - 



dVidVs \ftiri,r2)\ t^ifl,f2) 



(35) 



B. Dipole spin-spin interaction 

In the total spin representation, the dipole spin-spin inter- 
action can be rewritten as21 



Vdip — 



1.45 ^ / S'2r22 -3(S'-ri2 



P 2 



■ meV 



3 V^'-2)^(^-^ 



(30) 



Since S \xs) = and /t(fi = ra) = 0, where \xs) and ft are 
singlet-state spin and triplet-state coordinate wavefunctions. 



Taking into account the fact that the electrons are exponen- 
tially localized at sites rvi and tb in the /( state, a good ap- 
proximation to ti is to approximate the function ti by a con- 
stant value at those points where /t(ri, r2) is localized, thus 
obtaining the estimate 



TTtt 

-"dip 



0.36 



■ meV 



' AB 




(36) 



In order to further improve the estimate, the function ti (ri , r2) 
can be expanded in a Taylor's series around the localization 
points and the remaining integrals in the expansion terms be 
calculated analytically. From Eq. ( |36] | we find the estimate 
-ff^*p - (0.36/r^s)mcV w 5.0- 10"^ mcV at r^s = 100 i. 
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C. The _B-field interaction in the pure-spin model 

For the magnetic field one gets 



B, 



av 



+ 

av 



where 



Ba 



B„„ — 



B 

avx i iB, 



B- 

-Bn 



B 



(37) 



avyl 



(38) 



Using Eqs. ( fTTT i and (|37] |. one derives the Zeeman interac- 
tion Hamiltonian of the total spin S with the magnetic field 

Bat! • 



ij=i 



(39) 



Similarly, for the singlet-triplet matrix we have 

H^'iB) = 




~^SB- 
^SB+* 



(40) 



where 

AB 
SB 



B 



AB^ ± iAB„ 



B 



B 



SB^ ^ 6B,:±iSBy.(4-l) 



If the i3-field is homogeneous, from Eq. ( |4TI) we obtain 
AB = SB = and H'^B) = 0. In this case, the spin 
dynamics is unitary and is described by the Zeeman Hamilto- 
nian H**(i?) Eq. ( [39l ); the spin-spin dipole interaction iJ^-p is 
too small and can usually be safely neglected. 

Let us consider modifications due the to B-field inhomo- 
geneity in the /9Mre-spin model. Neglecting contributions from 
the double-occupancy states within the first-order perturbation 
approximation in the singlet-triplet interaction iJ"*, we find 
for the non-unitary term in Eq. ( fT4] i 



sin( J^fi) 
Jh 



([K,pKt 



[Kp,Kt]), 



where 



-^AB-J,, 



(42) 



(43) 



and 



S2 — Si X S 



= 2[S2 X Si]+Z(s2-Si) (44) 



S 0.7 
0.6 



Q- 0.5 
S 0.3 

Q. 

0.2 
0.1 



r =400 Angstrom 
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------0.30 

* 0.59 
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FIG. 1: (color online) The triplet states population of electrons in 
shallow QD centers in GaAs, as a function of time at different mag- 
netic field differences AB^ = 0.01,0.025, 0.05,0.075, 0.1,0. 2T, 
normalized to the Heisenberg exchange Jh constant. Initially the 
system is in the triplet state 15=1, Ms = 0). The distance between 
QDs is 400 A. 



is an asymmetric spin operator containing both linear and 
bilinear parts. Observe that Ct = at the "swap" times 
in — T^n/ Jh, n — Q,!, . . .. 

Similarly, for the Lamb shift in Eq. ( fT4l ) we have 



1 — cos( J/fi) 
1^ 



K^K. (45) 



Observe that Ct and Lt are quadratic in the difference field 
AB. Besides, notice that the magnetic field due to spin-orbit 
coupling does not contribute to the difference field Ai?so = 
but contributes to the i5i3so-field that is present in the coupling 
between the triplet states and the double occupancy, singlet 
states, in Eq. ( |40] i. If the external magnetic field Box is homo- 
geneous, then the singlet-triplet states coupling comes only 
from the spin-orbit interaction. Since the double-occupancy 
states should be involved in the dynamics in order to obtain 
non-zero spin-orbit interaction effects, these effects are ex- 
pected to be especially small, proportional to 5B^^, in the 
pMre-spin model. An estimate of these spin-orbit effects will 
be given in a numerical example in the next Section. 

Clearly, there is an important qualitative difference between 
pure- and pseudo-s^in models. In the former, the singlet- 
triplet states coupling is a second order effect, while in the 
latter this coupling is of first order in H^* [cf., Eq. ( l42l l and 
©J. Thus, in pMre-spin models effects due to B-field inho- 
mogeneity should be especially (quadratically) small as com- 
pared to the corresponding pseudo-s^m model effects. In case 
of negligible B-field inhomogeneity, as follows from Eq. ( l40b , 
the pMre-spin dynamics is unitary and is governed by the spin 
Hamiltonian H". 

Let us now consider a simple numerical example for the 
non-unitary effects due to the difference field AB for an elec- 
tron localized on a donor impurity in an n-doped GaAs semi- 
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FIG. 2: (color online) The purity p(t) = Tv p^{t) for the same 
parameters as in Fig. 1. 



FIG. 4: (color online) Triplet states population dependence on in- 
terdot separation rAB = 300, 350, 400, 450, 500, and 600 A, at a 
fixed AB^ = 0.05T. 
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FIG. 3: (color online) The Lamb shift energy as a function of time 
for the same parameters as in Fig. 1 . 



FIG. 5: (color online) Lamb energy shift dependence on interdot 
distance tab = 300, 350, 400, 450, 500, and 600A, at a fixed 
AB^ = 0.05T 



conductor. To simplify numerics, we assume that H'"' = 
disLg{esi,Ss2, ^ss) is diagonal and the singlet-triplet coupling 
field, AB, has only a non-zero z-component, AB^- Then, 
the corresponding eigenvalue problem, Eq. ( fT9] l. can be re- 
duced to a biquadratic polynomial equation which could, in 
principle, be solved exactly. If we neglect the exponentially 
small coupling field SB^, proportional to the overlap S, the 
biquadratic equation reduces to a quadratic one and we find 
for the non-unitary term 



_ and the Lamb shift 



K 



oj sin cut 



\ABl{l + cosujt) 
([K,pKt] + [Kp,Kt]) 

— -r^Bz ■ Jasz 



(46) 



Jh(1 — cos ujt) 
J| + iAB2(l + coswi) 



(47) 



where uj — ^JJj^+AB^. In the limit of small mag- 
netic field inhomogeneity, \ABz/Jh\ ^ 1, t^-" ^ \Jh\ 
and Eqs. (l46l ) and ( l47l ) go over into (l42l l and ( l45l l. respec- 
tively. Eq. WT\ describes the Lamb energy shift of the 
triplet state | S* = 1 , A/5 — 0) due to the coupling between 
singlet and triplet states induced by the magnetic field inho- 
mogeneity ABz- At the magnetic field geometry we have 
chosen there is no coupling between jS* = 1, A/5 = ±1) and 
IS- = l,A/s = 0) states. 

In Figs. [T][3] we show the results of calculations for the 
triplet states population, purity, and the Lamb shift energy, 
respectively, as a function of time at a fixed interdot separa- 
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tion (rAB = 400 A) and different AB^- For the Heisenberg 
interaction constant we used an asymptotically correct 
expression22i22i^ obtained for hydrogenlike centers in GaAs 
[note that our Jh — et — Eg is related to the exchange in- 
tegral J in Ref. 28 via Jh = — 2J]. Initially, the system 
is assumed to be in the |5 = 1, Ms — 0) state. As can be 
seen from Fig. 1, there is a re-distribution between singlet 
and triplet states population due to the singlet-triplet sub- 
space coupling. At \ABz/Jh\ ^ 0.1, the probability of re- 
distribution is negligible and the time-evolution is basically 
unitary. With increasing AB^, this probability re-distribution 
is seen to be more pronounced, time-evolution becomes non- 
unitary (Fig. 2), and |at(t)p can drop to the value Jjj/oj'^ at 
t = TTn/uj, n = 1,3,.... Observe that the non-unitary dy- 
namics reveals repetitions in time and at moments of maximal 
(minimal) singlet-triplet states probability re-distribution we 
find maximal (minimal) Lamb energy shifts (Fig. 3). Thus, 
the non-unitary effects observed are not irreversible and they 
do not result in a real decoherence process. We do not have 
in our two-electron model a real, external and infinite "bath", 
coupling to which would result in irreversible decoherence ef- 
fects in the spin system. In Figs. 4,5 we demonstrate the de- 
pendence of triplet states population and Lamb energy shifts 
on the interdot distance tab at a fixed AB^ = 0.05T. 



D. Spin-orbit interaction in pure-spin model 

In this subsection we estimate the non-unitary effects in 
the pure-spin model due to spin-orbit interaction. For sim- 
plicity we assume that the external magnetic field is homo- 
geneous and directed along the z-axis, with Boz being its z- 
component. Since 6Bso is a pure imaginary field (its compo- 
nents are matrix elements between the real states (pA and (j>B 
of an odd vector function of the momentum operator, both in 
vacuum and in the bulk of semiconductors that lack inversion 
symmetry, Dresselhaus fields^ as well as in heterostructure 
zinc -blendes, Rashba fields, ^^) we have 



-SBfo, 



6b: 



-SB, 



(48) 



Using these relationships, the singlet-triplet spin-orbit cou- 
pUng can be written as 















H^'iBso) - 



lSB+^ ~-^6Bsoz -\6B'^^ 



(49) 



The couplings between double-occupancy, singlet and 
triplet states are seen to be the same. We assume that 
i?'"* = dia.g{£s,edo,£do), where and e^o are the sin- 
glet and double-occupancy states energies, and — 
diag(et_|_, £t-), where £t is a triplet state energy and £t± — 
£t ± Bqz- Within these approximations, the 6-by-6 eigenvalue 
problem Eq. (fT9] l is then reduced to computing the roots of 
the biquadratic equation^'"' 



asE^ + 02 £^ + aiE + oq = 



(50) 



where 



03 



4 
i=l 



02 



E m 



2 

soa I 5 



a—x,y,z 



\5Bs 



\[\5Bs 



\5Bsoy\''] ) (£2+64) 



£i 



{\5Bsox? + \5Bsoy\^)£3, 
4 

= - |(5Bso2|^£2£4 - 



1=1 

\[\5Bso.\ 
£s, £2 — £t 



\5Bsoy\ ]£3(£2 + £4), 

£3 = £*, 



£4 = £t- 



For hydrogenlike centers one can estimate the energies e^o 
and £t as follows. The ground energy of two well separated 
hydrogen atoms is i?2H ~ — 27.2eV . Using for GaAs the 
scaling factor i^caAs = m*/£^ « 4.6 • 10^^ one can estimate 
£t ~ KqiiAsE2h — — 12.6meV. £do is located higher than 
£t due to mainly interelectron repulsion Vco{AA; AA) so that 
£do -£t = Vee{AA; AA) « 12.6 mcV. 

If SBsoa = 0, a — x,y, z, the roots of Eq. ( l50l ) Ei are 
equal to £i, i ~ 1, . . . , 4. The two other roots are E5 — £do 
and Eq = £5. The corresponding eigenvectors are 



= d, 0,1,1,- 



SB- 



V2SBs 



SB^ 



E, 



£2 ' Ek ~ £3 ' Ek — £4 



k = 1,...,4, 

leg) = (0, 1, -1,0,0, of , 

lee) = (1,0,0,0,0,0)^, 
where 



(51) 



dk = 2 



[\SB,o.\ 
1 



\SBsoy\ 
1 



iEk-£2y {Ek-£if 



2\5B, 



{Ek - £3? 



^1/2 



Notice that the above formulas are not valid in the degener- 
ate case: Bqz = and £2 = £3 = £4 = £t- In this case the 
biquadratic Eq. ( fSOl l reduces to two quadratic ones, two roots 
of which are degenerate, Ei — E2 — £t- Formally, one gets 
singularities in Eq. ( ISTI l at Ei = E2 = £t- Therefore, the 
simpler, degenerate case should be analyzed separately and 
the corresponding formulas [not shown here] can be derived. 

Let us now find the spin-orbit field 



SBso = [4>A B,M ^b)~ {MBso{p)\c^b) 

= j dr4>A{\r-R\)Bso{-i'^r)(t>B{r) (52) 
= Bso{-iV^)Sdr4>A{\r-R\)Mr) 
= B,,i-iVj^)SiR) 
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FIG. 6: (color online) Triplet states population dependence on spin- 
orbit interaction as a function of time, at different interdot separations 
Tab and fixed orientation of the interdot radius-vector, tab, {d = 
7r/4, 1^ = 7r/3) [see text]. The initial triplet states population is 
taken to be equal. The external magnetic field Boz = 0.5T. 



FIG. 7: (color online) The purity in the presence of spin-orbit cou- 
pling. All parameters are the same as in Fig. 6. 



where Bso{p) is an odd function of the momentum operator 
p — — iVp and R = tab- In particular, in zinc-blende semi- 
conductors such as GaAs, Bgo is cubic in the components of 

p: 



'■31.34 



{cyclic permutations of x,y,z\ 



Bsoa 



'2m*E, 



(53) 



where m* is the effective mass of the electron. Eg is the band 
gap [m* w 0.072, Eg k, 1.43 eV for GaAs], px, Py, Pz are 
components of the momentum along the cubic axes [100], 
[010], and [001] respectively. The dimensionless coefficient 
aso = 0.07 for GaAs. From Eqs. ( l52b and ( |53] ) we obtain 



Bs 



iAs 



i?3 



{s"'iRy 



y\R) + l^S'{R) 



The overlap integral-'^ S{r = R/as) = (1 



(54) 



r /3)exp(— r) for hydrogenlike centers [as ~ 92 A for 
GaAs] and Eq. (l54l i reduces to 



Bsox — i(0.83meV) sin f? COS (^(sin^ 6*8111^ — cos^ 



(55) 



where {R, 6, ip) are spherical coordinates of the vector R, with 
other components being obtained by cyclic interchange of in- 
dices. 

In Figs. |6] and |7] we display the time-dependence of the 
triplet states population and the purity, which is induced by 
the spin-orbit interaction, Eq. ( fSSl l, at a fixed orientation. 
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FIG. 8: (color online) Effect of spin-orbit angular anisotropy on 
triplet states population. All parameters are the same as in Fig. 6. 



(9 = 7r/4, </? = 7r/3), and different r^B in the range 200- 
500A. Observe that the maximal re-distribution of singlet- 
triplet probability occurs at 27rt/|erfo—et| = 0.5 and the spin- 
orbit interaction effect diminishes as tab increases. The max- 
imal singlet-state probability achieved at tab — 200 A is seen 
to be quite small, ^ 10^''. As compared to the non-unitary ef- 
fects induced by S -field inhomogeneity, the spin-orbit effects 
are on average four orders of magnitude smaller. The angular 
dependence of the population of triplet states on the interdot 
radius-vector orientation at a fixed tab — 200 A is illustrated 
in Fig. [8] 
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E. The _B-fleld interaction in the pseudo-spin model 

Using Eqs. JJTb and ( l40l l. the effective Hamiltonian matrix 
Eq. (|6|l in the basis {$sii ^u, * = 1, 2, 3} can be rewritten as 



Ti-eff — 



2^2 



2V2 2 ^ 2\/2 







av 

B+ 



^ J 

(56) 



where for simplicity we neglected contributions from the 
double-occupancy states [the resolvent term in Eq. ([3])]. Al- 
ternatively, in the pseudo-spin representation we get 



1 



'HcB — I -£ 



B„ 



, + -Ss j I + Jh Sa ■ SB + 

- N ^B 

( SA + Sb) + — - (SA- Sb) 



1 3 

— F -\ 

Ba ■ SA + Bb ■ SA, 



St] I + Jh SA ■ SB 



(57) 



where Ba = Bav + AB/2 and Bb = Bav - AB/2 are the 
local magnetic fields at sites A and B, respectively. The term 
Jh Sa ■ SB is the familiar Heisenberg interaction. In matrix 
form, Eq. (|57] | can be rewritten as 



"Wcff — 

ei 
£3 



Bl 


V 

1 

2 
1 



St 



\Ba 

£2 

hJn 



Baz 





\Jh 

1 R^ 



2^B 

£2 = £t 






£4 



1 



/ 

( + Baz 
= £. - 



(58) 



where B 



A,B 



Ba,Bx ± iBA,By 



The Hamiltonian Eq. ( |56] ) generates a unitary evolution 

t/eff (t) = exp(-imeff) (59) 

in C*. At a fixed set of parameters £s, £(, Ba, Bb the propa- 
gator Ucs{t) does not provide a universal set of unitary gates 
in C"'. Any unitary transformation U G U{4:) can be repre- 
sented as a product of a phase factor exp(iQ;), where a is a 
real parameter, and a unitary transformation Us G SU{4). 
Any transformation t/5 is determined by M = 4^ — 1 = 15 
independent real parameters (6'i, . . . , 615) so that 



UsiOi,. 



) = exp 




(60) 



where the set of generators {Fi] is an orthonormalized trace- 
less, Hermitian matrix set that forms a Lie algebra su(4) [F^'s 
form a complete basis in a real M-dimensional vector space; 
they are analogs of Pauli matrices, ctq, a = x,y, z, in su{2) - 



see, e.g., Ref.'s^. From the representation Eq. ( l60b it follows 
immediately that UcS,s{t) = exp (-it (ie^ + |ef)) Ucs{t) 
cannot match an arbitrary Us, because the number of inde- 
pendent parameters in Eq. (l56b is at most 8 - fewer than the 
number of di's. This can also be understood from the fact that 
the form of the Hamiltonian matrix, Eq. ( |56] |, is not generic. 
In particular, the matrix is sparse, i.e., the entries (2,4) and 
(4, 2) are zeros. 

However, compositions of unitary transformations Eq. (|59] l 
taken at different sets of parameters can provide a universal 
set of unitary gates in C"'. A well-known example of universal 
gates is provided by the Heisenberg interaction (at Jh 7^ 0) 
with single-spin addressing (at AB 7^ 0).' 

From Eq. (|56] l it follows that a necessary and sufficient 
condition for obtaining a universal set of gates on two spins 
is to have an inhomogeneity in the magnetic field AB ^ 
[the source of inhomogeneity can be different, it can be ei- 
ther strongly localized magnetic fields or g-factor engineer- 
ing] and the Heisenberg interaction, Jh 7^ 0. The reason 
is that when AB = 0, the Hamiltonian Eq. ( |56] | and the 
corresponding unitary transformations take a block-diagonal 
form, with singlet-triplet entries being zeros, while when 
Jh = 0, the Hamiltonian form (l58l l will have zero off- 
diagonal block-matrices. Clearly, even a composition of such 
unitary transformations taken at different sets of parameters, 

either £s,£t,BA — Bb or Es — St, Ba, Bb, will be in a 
block-diagonal form and it cannot reproduce an arbitrary uni- 
tary transformation. Note that when one allows for encoding 
a qubit into three or more spins, the Heisenberg interaction 
alone is universal in the pseudo-spin model,^^'^^ and Heisen- 
berg along with an inhomogneous magnetic field is universal 
for an encoding of a single qubit into pair of spinsi^ 

Moreover, it should be noted that in the homogeneous mag- 
netic field case unitary transformations restricted to the triplet 
subspace will not provide a universal set of gates. To prove 
this statement, let us consider a composition of two unitary 
transformations in the triplet subspace: 



exp(-iiiF"(Bi)) exp(-it2ff"02)) = 
exp(-i(ti +t2)ifcff) = 

(-t{ti + t2)£t - itiH'\Bi) - it2H'\B2y 



exp 

Ma 
2 



(61) 



H'\Bi),H'\B2 



where on the right-hand-side we used the Campbell-Hausdorff 
formula.'*'^ From Eq. ( |37] ) one obtains 



H'HBi),H'\B2) 




(62) 



Bx ± iBy. Since the 



where B =[Bi x B2] and B^ 
higher-order terms in the Campbell-Hausdorff formula con- 
sist of nested commutators between H*'^{Bi) and H^^{B2), 
we find that the effective Hamiltonian 77*^ corresponding to 
the product of two unitary transformations will still have a 
sparse form, with the (1, 3) and (3, 1) entries being zeros. 
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F. Is it possible to obtain a universal set of gates in the 
pure-spin model? 

Simultaneously, we have just proved that in the pMre-spin 
model, in the case of a homogeneous magnetic field, unitary 
transformations in the triplet subspace will not provide a uni- 
versal set of gates. On the other hand, at AB ^ 0, we have 
already shown that the evolution of the spin-density matrix 
is non-unitary. Let us assume that we have a non-unitary 
gate £i so that p{ti) = Ci{ti)p{0). How could one define 
a composition of two non-unitary gates, C2C1I In order to 
do this unambiguously, £2 should obey a compatibility con- 
dition with the initial state [because a non-unitary £-gate is 
not totally independent of the initial state, it includes some 
sort of correlation information encoded in the initial state], 
that is a correlation Rm{ti) established between as(ti) and 
at{ti) amplitudes ait = ti should be included in in the def- 
inition of the corresponding dynamics generator operators in 
£2- Eq. (fTTI i, where the left-hand-side and H*^ should be 
replaced by R„i{ti) and the identity matrix, respectively, pro- 
vides a relationship between Rm{ti) and i?m(0). If the corre- 
lation between the amplitudes at i = and t — tiis the same, 
Rm{ti) — Rm{0), then we obviously have Ci = £2 — ^ and 

£2it2)Clitl)=C{ti+t2). 

In the total Hilbert space, the state is defined by 1 1 real pa- 
rameters. While in the reduced description, the spin-density 
matrix is defined by 5 real parameters (for more on the spin 
density matrix parametrization in terms of a's, see SectionlVb. 
Fixing a correlation in the initial state, as{0) — i?m(0)a((0), 
we have 3 complex equations between the amplitudes 0^(0) 
and af(0), which define a 5D real manifold embedded into 
the total Hilbert space. Using these equations we can separate 
6 extra real degrees of freedom that we have in the total state 
description from those in the spin-density matrix description. 
However, these extra degrees of freedom are not eliminated 
in the spin-density description, they are included in the form 
of correlation matrix Ra{0), a = {s, t, m}. It was shown in 
part I that Eq. (fT4] l provides an exact description of quantum 
evolution, in the spin-density space. Therefore, as long as we 
have a universal set of unitary gates in the total Hilbert space, 
this set of gates will be translated into the corresponding uni- 
versal set of non-unitary gates generated by Eq. ( fT4l l because 
no information is lost in our "reduced" spin-density matrix 
description. 



G. Spin-orbit interaction in the pseudo-spin model 



Let us consider spin-orbit effects, which are proportional to 
6Bso, in the pseudo-spin model. From Eq. (|6]l we obtain 



where 



kk'=2,3 



1 



EI- H'' 



Fk'iSBso), 



kk' 



FkiSB^o) = J2 (Hh{SB,o)K, + {H'^{6B,,)rKl 



(64) 

It follows from Eq. ^ that = i?!*, F2 = F3 and 
Eq. ( l64l l can be reduced to 

F2,3{^Bso) = F{6Bso) = iV2SB,o [sa x sb] (65) 
Then, Eq. ( |63] l can be rewritten as 



nsoiSB,o)^AdoF{SBso) 



(66) 



where the coefficient Ado is proportional to the ratios of the 
amplitudes of double occupancy transitions (i?i| and Hf^) 
and the energies of interelectron interaction (iJ|| — E and 
— E) in doubly occupied QDs: 



Ad 



kk'=2,3 

Hl^{E^ 

TTSS 

E - H^. 



1 



EI- H'^ 



kk' 



12 (F - 



^32)+ 



m 



E- H: 



33 



Here, the determinant is 



A^{E-H^I)iE-H^I)-\H, 



ss \2 



(67) 



(68) 



The effective spin-orbit interaction Hamiltonian Eq. ( I66b is 
different from the corresponding one obtained by Kavokin.'- 
In our derivation the double-occupancy states are essential, 
whereas in Ref. 17 these states are totally neglected. We 
showed above that neglecting double-occupancy states results 
in zero spin-orbit coupling. The very physical picture put for- 
ward in Ref. 17 to support the derivation was based on the as- 
sumption that when one of the two electrons localized at cen- 
ters AoT B tunnels to the adjacent center (say, from A to B), it 
experiences the influence of the spin-orbit field resulting from 
the under-barrier motion of the electron. Neglecting double 
occupancy states means that the second electron should si- 
multaneously tunnel from i? to A so that the two electrons can 
never be found in the same QD. Indeed, in Ref. J/7. this simul- 
taneous two-electron transition is described by the product of 
two matrix elements: the overlap S and {SBso)a, ce = x,y,z 
[Hso ^ S{6Bso ■ [sA X Sb] ), our SBgo is related to Kavokin's 
b-field via 6Bso = —ih, and the overlap via 5 = il]. With 
the orthogonalized molecular-type, two-center orbitals such 
a one step two-electron transition gives a zero contribution 
since the spin-orbit interaction is a one-electron operator and 

the overlap S = 



(63) 



'JB ^ — 0. Eq. i63[ describes the 

two-step mechanism: in the first step, the two-electron sys- 
tem makes a transition from the singly-occupied state $^1 to 
the intermediate, double-occupancy states $32, *f s3 due to the 
inter-electron interaction (Hf^, k — 1,2 terms). Then in the 
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which is exactly the ratio of our and Kavokin's estimates, 
as a function of interdot distance. As one can see, Kso de- 
creases from 0.98 to 0.46 in the range tab ^ 300 - 700A. 
Interestingly, our results qualitatively agree with the results 
of Ref. |2^ which obtained in the region of interest [tab ~ 
(3 — Tlanl a reduction of about one-half relative to that of 
Refs. T7'!4r. According to Ref. 28, Hso ^ 2 J7GK, where J 
is the exchange integral calculated using the medium hyper- 
plane method22i2£ and 7gk is an angle of spin rotation due 
to spin-orbit interaction introduced by Gor'kov and Krotkov 
[7gk ~ hl'K^ 7k being the corresponding angle of spin rota- 
tion introduced by Kavokin]. Note that Kavokin's 7k and J 
are not independent parameters. In Ref. XL 7k was defined 
as Qb/ J so that their product 2 J7K = 2ri& does not depend 
explicitly on J. 



FIG. 9: Spin-orbit interaction reduction coefficient Eq. J73t as a 
function of interdot distance for GaAs. 



IV. ELECTRONS ON LIQUID HELIUM 



second step, as a result of the spin-orbit interaction, the sys- 
tem makes transitions from $32, to triplet states (the 
Hll terms). 



Let us find an estimate for 
2V2 



Ado{R) 



\£do — £t\ 



- R\) 



dridr2 (f>'^{ri) x 



1 



e\r2 - ri 



■0(^2), 



(69) 



where the hydrogenlike orbital (j){r) 
(ttc'^)^^/^ exp(— r/as). Since electron 1 in Eq. 
is localized around the effective Bohr radius qb, one can 
approximate 



1 



1 



\r2 - ri\ 



\f2 - aBn/nl' 



(70) 



Then, the remaining integrals can be calculated exactly and 
we obtain 



Ado{r = R/aB) 



-i^{4eFi(r)- 
2[F,{r + l) + F,{r-l)]- 
[F2{r + 1)-F2{r-1)] + 
2 [Fo{r + 1) - Fo{r - 1)]} 



(71) 



signjx) 



[I5{l + |x|) + 6|x|2 + |x|3] exp(-|:E|), 



where 

Foix) 

F,{x) = ^[3|x|(l + |x|) + |x|3]exp(-|x|), 
F2ix) = SiS£(£)(i + |a;|)exp(-|x|). 

(72) 

Here, we used the same estimate for £do — £t as in Section 

mm 

In order to compare our calculations to Kavokin's result for 
GaAs, we have plotted in Fig. |9] the spin-orbit interaction re- 
duction coefficient 



\Ado\ 
V2S' 



(73) 



Recently, Lyon suggested that the spin of electrons floating 
on the surface of liquid helium (LHe) will make an excellent 
qubit.**^ Lyon's proposal, instead of using the spatial part of 
the electron wavefunction as a qubit as in the charge-based 
proposal ,^^'"^'"^^'*^^ takes advantage of the smaller vulnerability 
of the electron's spin to external magnetic perturbations and, 
as a consequence, a longer spin-coherence time. It also has the 
important advantage over semiconductor spin-based propos- 
als (as first pointed out in the charge-based proposal'*-''^'*''*'''^-) 
that, with the electrons residing in the vacuum, several im- 
portant sources of spin decoherence are eliminated so that the 
environment effects are highly suppressed (the spin-coherence 
time is estimated^ to be T2 > 100s). 

The geometry of the system with the electrons trapped at 
the LHe-vacuum interface (see Fig. 2 in Ref. 42), is concep- 
tually similar to that of semiconductor heterostructures. The 
two electrons are trapped at sites ta and re {rA,B are radius- 
vectors of the centers of quantum dots in the z = plane) 
by the two attractive centers created by two charged spherical 
electrodes located below the LHe surface at distance h and 
separated by the interdot distance tab 



rA - rB 



A 



Vtriz,7^ = -- + VA{z,r} + VB{z,f} (74) 



VA,B{z,f) = - 



Qa,i 



{z + hf + {r-fA,Bf 



(75) 



where the first potential, — A/z, is due to attraction to the im- 
age charge induced by an electron in the LHe [A = (e — 
l)/(4(£ + 1)) w 7 • 10"^, with e w 1.057 being the dielectric 
constant of helium] . For purposes of interaction control, the 
Qa.b charges on the electrodes can be made variable in time. 
The electrons are prevented from penetrating into helium by a 
high potential barrier (^ 1 eV) at the helium surface, so that 
formally one can put Vtr = 00 at z < 0. The in-plane and 
out-of -plane motion of electrons in the potential of Eq. ( |74l i is 
in general non-separable. However, near the electrode's posi- 
tion 'rA,B the potential of Eq. ( iTST i is approximately separable 
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in the z and r coordinates 

I- -^ ^ Qa,b , 1 2 - ^2 

Vk,S f'j ~ V C^A.B Z + 2^A,B {l~ - ^A,B j , 

(76) 

where it is assumed that z and \r — rvi.s| ^ h and = 

1/2 

Q A,B I h"^ , ^ A,B — {QA,B/h^) ■ In the separable approx- 
imation of Eq. (l76l l. the electron's motion in the z-direction 
[z > 0] is described by a ID-Coulomb potential perturbed 
by a small Stark interaction and the in-plane motion by a 2D- 
oscillatory potential. We assume that the out-of-plane mo- 
tion in the z-direction is "frozen" in the ground state of a ID- 
Coulomb potential 

ip^{z) = 2VA(Az) exp(-Az), (77) 

and the in-plane motion, in a superposition of the in-plane 
confining oscillatory potential and, possibly, a perpendicular 
magnetic field, is described by the Fock-Darwin (FD) states 
of Eq. ( l23b . Then, the calculation of hij and Veeiij'jkl) in 
Eqs. (|27]i-(|29l) in the chosen basis set is reduced to the one- 
dimensional integrals 




500 600 700 800 900 1000 

Interdot Distance r^^, Angstrom 



FIG. 10: The Heisenberg and dipole-dipole spin interaction mag- 
nitudes as a function of interdot distance. The distance from the 
electrodes to the LHe surface h = 800 A and the charges on the 
electrodes Qa = Qb = 1 a.u. The magnetic field -Bo ~ 0. 



h-Ai 



hAA 



T-BB 



A 
2" 



-2AQAgc{a,(3,XA)- 



2A(3_Bffc(a, /?, As)] Sab, 
-| + ^ - 2AQAgc{a, Pa, 0) - 
2AQBgcia,(3A,'2ArAB), 
- - ^ ^ - 2AQAgcia,(3B,2ArAB) 



41 



B 



5c(a,/3, A) 



2AQBgc{a,pB,0), 

/"OC 

/ dx Jo (Ax) X 
Jq 

exp {-ax - I3x^) /{x + 1)^, 

{2AlAlBf 



(78) 



= 2Ah, /3 = 



/2 



/2 



PA,B = 2[AIamYi ^a,b 



2Al%, 



Y'^AB, 



Vee{AB-CD) 



gee{a,b) 



Neegeei.a,b), 

AIaIbIcId 

^2 

exp ' ^'^ 



' BD 



4{l\ + ll) 4(1% + 11) 

dx Jq (bx) X 
(3x2 + 9x + 8)(x + l)"^ 



(79) 



= 4A^ 



2A 



72 /2 



72 72 

''B'-D 

72 I 72 



Ic^A + I'aTc l^rB + l^ro 



72 I 72 

'B ^ '■D 



where i = A,B,C,D in the two-electron matrix elements 
denotes orbitals with the effective lengths U, localized at r^, 
and Jq{x) is the zeroth order Bessel function. From Eq. ( |79l ) 
one can obtain the following expression for the Heisenberg 
interaction constant 



Jh 



-1.36- 10^A5ee(a,0)S'2meV 



(80) 



72 I 72 ■ 



Note that Ju is proportional to the square of the overlap ma- 
trix element S, Eq. (IZST i, and the integral gee{a, 0) does not 
depend on the interdot distance tab- As a rough estimate, one 
can approximate the rational function in the integral gee{0', 0) 
by a constant 8, as a result obtaining gee (a, 0) « Ay^Tr/a. 

Figure [TO] shows the magnitudes of the Heisenberg and 
dipole-dipole spin interaction, \Jh\ and Jdip from Eqs. dSOl ) 
and ( |36] | respectively, as a function of interdot distance. As 
estimated, the magnitude of the Heisenberg interaction is 
comparable to the weak dipole-dipole interaction at tab — 
1000 A. However, we remark that the strong dependence of 
S ~ exp(— ar^^) (quadratic in the interdot distance) is due 
to the quadratic dependence on coordinates in the exponent of 
the corresponding oscillatory wavefunctions. Asymptotically, 
the confining potential Eq. dTSb behaves as a 2D Coulomb 
potential, so that one should expect a milder coordinate de- 
pendence, S ^ exp(— ar^s), at large distances (assuming 
Bq = 0) and the rough estimate of Eq. dSOb provides a lower 
bound for the Heisenberg interaction strength. 

Lyon suggested"*^ using, instead of the exchange interac- 
tion, the magnetic dipole-dipole interaction between the spins 
in order to implement two-qubit gates, motivating this by the 
strong sensitivity of the exchange coupling to the parameters 
of the system and, hence, the corresponding difficulties with 
attempting to control this interaction. Our analysis confirms 
this, though, of course it is not easy to control the dipole- 
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— I ' 1 ' 1 ' 1 ' 1 ' r- 

0.0 0.2 0.4 0.6 0.8 1.0 

time, 2jtt/|J^| 



FIG. 11: (color online) The triplet states population of electrons 
trapped in QDs on a LHe surface, as a function of time at differ- 
ent magnetic field differences ABz = 0.05 and 0.1 T. Initially the 
system is in the triplet state |S = 1, Ms ~ 0). The distance between 
QDs is 900 A. All other parameters are the same as in Fig.llOl 




— 1 ' 1 ' 1 ' 1 ' 1 ' r— 

0.0 0.2 0.4 0.6 0.8 1.0 

time, 2kV\JJ 



FIG. 12: (color online) The purity p{t) for the same parameters as in 
Fig.E] 



dipole interaction either: Eq. ( |36] | shows that this interaction 
depends on only one controllable parameter, the interdot dis- 
tance (the 5-factor is a constant in vacuum). 

Similarly to Figs. [T] and |2l Figs. [TT] and [12] demonstrate 
non-unitary effects in the pure-spin model, due to magnetic 
field inhomogeneity in the electrons-on-LHe system. The in- 
terdot distance shown is rAB = 900 A. At this distance the 
Heisenberg interaction still prevails over the dipole interac- 
tion by at least an order of magnitude. Again, the pattern seen 
in the singlet-triplet states population redistribution is clearly 
oscillatory. 



V. DISCUSSION AND CONCLUSION 

We have performed a comparative study of pure- and 
pseudo-spin dynamics for a system of two interacting elec- 
trons trapped in two QDs. We have shown that when there 
is negligible coupling between the spin and orbital degrees 
of freedom, which is the case of near i?-field homogeneity 
and negligible spin-orbit interaction, the system spin dynam- 
ics is unitary in both pure- and pseudo-spin models and is 
governed by the Zeeman interaction Hamiltonian of the to- 
tal spin S (S = 1) with the magnetic field Bav The sin- 
glet and triplet states are totally decoupled; the total spin is 
conserved. The spin system Hilbert space can be decom- 
posed into two independent, singlet and triplet subspaces, 
the singlet spin states being magnetically inactive (5* = 0). 
Thus, the two-electron spin system restricted to the triplet 
subspace physically embodies a qutrit. The Heisenberg in- 
teraction operates differently in pure- and pseudo-spin mod- 
els. If for simplicity we neglect double-occupancy states, the 
pseudo-spin state is totally defined by four complex ampli- 
tudes: {asi{0),au{0), i = 1,2,3} in the basis {^si,^u}; 
so that the Heisenberg interaction results in a phase uni- 
tary transformation: {asi{t) — exp{—iest)asi{0), ati{t) = 
exp{~iett)au{Q)}- Since the spin-density matrix p{t) is a 
bilinear combination of a's: ptij{t) — ati{t)a1^{t), Ps{t) — 
1 — Ptii(t), the yO-state will not by affected by the unitary 
transformation induced by the Heisenberg interaction. 

We have also shown that unitary quantum gates realized in 
both spin models do not provide a universal set of gates un- 
der the condition /S.B = 0. In order to obtain a universal 
set of gates, there should be both non-zero coupling between 
singlet and triplet states (AB ^ 0) and non-zero Heisenberg 
interaction ( Ja ^ 0). Although at /S.B ^ pure-spin dy- 
namics becomes non-unitary, one can establish a relationship 
between unitary gates in pseudo-spin and the corresponding 
non-unitary gates in pMre-spin dynamics so that a universal set 
of quantum gates constructed within the pseudo-spin model 
will generate a universal set of non-unitary gates in pure-spin 
dynamics. 

To demonstrate the non-unitary effects, which are pro- 
portional to the square of the magnetic field inhomogene- 
ity, in pure-spin dynamics, we have calculated how singlet 
and triplet states populations, as well as the purity and the 
Lamb energy shift, are affected by AS ^ and spin- 
orbit interaction in n-doped GaAs semiconductors. These ef- 
fects are found to be strongly dependent on the ratio of B- 
field inhomogeneity and the Heisenberg interaction constant, 
\lS.Bz/ Jh\- For example, the singlet-triplet states popula- 
tion re-distribution is maximal at t = vrn/w, where w — 
^/j^H-ASf and n = 1,3,..., and the singlet state pop- 
ulation can achieve the value /S.B'l/uP'. Thus, we can con- 
clude that the Heisenberg interaction, characterized by the 
interaction constant Jh, plays an essential role in producing 
non-unitary effects in pure-spin dynamics. Spin-orbit interac- 
tion effects are found to be roughly four orders of magnitude 
smaller as compared to those caused by B-field inhomogene- 
ity. 
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As shown in Figs. [T][8] there are clear oscillations in the 
pure-spin dynamics and the non-unitary behavior of the spin- 
density matrix does not show the decaying pattern character- 
istic of a real decoherence process. This should be expected 
since the "bath" - the electron orbitals - in our spin model 
is not a real stochastic or infinite external bath, an interaction 
with which may result in irreversible decoherence. In essence, 
the spin dynamics is embedded in space and our bath is too 
small. The coordinate Hilbert space in the two-orbital ground 
state approximation adopted in the present paper is repre- 
sented by 4 two-electron coordinate basis wavefunctions. In 
principle, the coordinate bath can be large in a system where 
couplings between excited and ground state orbitals are not 
negligible. This is an interesting question for future investi- 
gation: how will couplings to excited orbitals affect the non- 
unitary spin dynamics? The other interesting generalization of 
the present model is inclusion of real environment effects, i.e., 
the real stochastic bath representing the interaction of electron 
spins with the semiconductor medium. We will consider these 
and other generalizations in a future publications. 

In the pseudo-spin model, where B-field inhomogeneity 
results in first-order effects, we have estimated the contri- 
bution of the spin-orbit interaction to the effective pseudo- 
spin Hamiltonian, namely the Dzyaloshinskii-Moriya spin- 
orbit interaction term, and have suggested a two-step mecha- 
nism: coupling between the singly-occupied singlet state and 
triplet states occurs via intermediate, double-occupancy states 
(direct coupling between these states turns out to be zero due 
to orthogonahty of the orbitals involved in the transition). Our 



calculations predict a smaller magnitude of the spin-orbit in- 
teraction as compared to the estimates of Ref. 17, but are con- 
sistent with the results of Ref.l28i 

In our second application we demonstrated, in Figs. [TTI and 
[T2I non-unitary effects due to AB 7^ in a system of elec- 
trons trapped above a liquid helium surface, namely the spin- 
based quantum computing proposal by Lyon^^. A more thor- 
ough investigation of spin dynamics in this system is left for a 
future publication. 

In conclusion, we note that the two-electron spin-density 
matrix description advocated in this paper is expected to be 
useful when electrons trapped in QDs are not spatially re- 
solved or resolvable. The spin-dynamics is then completely 
described by the spin-density matrix p. Although the p- 
dynamics becomes non-unitary in general (at AB ^ 0), it 
is controllable by modulating the interaction parameters, Jh, 
Bav, AB. Since the non-unitarity comes from the mag- 
netic field inhomogeneities, AB and/or 6Bso, and since the 
p-dynamics is quadratically protected from these fields, this 
might prove to be important in practical quantum computing 
as minimizing coupling between spin and orbital degrees of 
freedom quadratically improves the fidelity of unitary gates 
in the p-state space. 
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